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Learning Aims: You will learn to 



State and use the phonon energy distribution. Derive the 
density of states for phonons. Explain and derive the Debye 
frequency. 

Derive the low and high temperature limits of the phonon heat 
capacity. 

Explain blackbody radiation as a free photon gas. State and 
use the photon energy distribution. Derive the density of states 
for photons. 

Derive Stefan's law for blackbody radiation. 



Specific heat capacities of solids 



In 1912, Peter Debye, a Dutch physicist working in Germany, 
produced a theory which predicts the heat capacity of a solid 
correctly at high and low temperatures. 
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Atoms vibrating in a solid 



In Debye's model, we think of the atoms in a solid as 
connected by springs. They vibrate in 3D in a complicated way. 

The problem can be simplified by treating the vibration as a 
superposition of waves of different frequencies. 




If we suppose that the atoms are fixed at the edge of the solid, 
the frequency would be quantised. A quantum of this vibration 
energy is called a phonon. 
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Phonon occupation number 



The phonon energy distribution is given by 

n{e)de = ; — , r 

exp(e/fe B T) - 1 

Next, we need a formula for the density of states g{e). 

Note that although we have previously used the same symbol 
g{e) for the particle in a 3-D box, the phonon density of states 
would have a different formula. 
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Quantised thermal waves 



We can count the number of states for the phonons, like we 
did for atoms in an ideal gas, or electrons in a metal. 

Unlike electrons, however, the exclusion principle does not 
apply to phonons. 

So each energy state can be occupied by any number of 
phonons. 



density of states for phonons 

The vibrations in a solid is described by a wave equation. The 
solution for the displacement of the atoms £ , can be given by 
the same function 

£ = sin k x x 

in each of the x, y and z directions. Assuming that the 

displacements are zero at the edge of the solid, the 

wavevectors are discretised in the same way: 

_n x iv 
n x — ■ 
a 

So the density of states can also be given by the same formula: 

Vk 2 

However, two changes are needed from the way we have used it 
for the ideal gas and the electrons. 



#)- 
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density of states for phonons 



The first change is in the change of variable to energy e. 

For atoms or electrons, the energy is 

h 2 k 2 



For vibrations, the energy is given by by the same formula as 
for photons: 

e = Tkj, 

where uj is the angular frequency 2nf and / is the frequency, uj 
is in turn related to k by 

UJ 

where v is the velocity. 



About polarisation 



The second change is the related to the polarisation. Recall 
that for electrons, we have to multiply by 2 because of the 2 
spin states. 

For vibration, a state is a particular frequency, for a particular 
set of wavevectors (k x , k y , k z ). 

Each state correspond to a wave, and a wave in a solid can 
have 2 transverse polarisations, and 1 longitudinal polarisation. 

So we have to multiply by 3 to get the correct density of state: 

VuJ 2 duj 

g{uj)duj = 3 x 
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density of states for phonons 



The variable of angular frequency w is often used instead of 
energy e. For the change of variable, we then use 

oj = vk. 

We can now do the change of variable using 

g(uj)duj = g(k)dk. 

Substituting 



fl(*0 - 



Vk 2 
2tt 

we get the density of states in uj: 

< s Vu}2 

9 ^-> = o 2 3 ' 
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The Debye Frequency 



There is one more practical point to note. In a solid, there is a 
limit to the highest phonon frequency. 

This is related to the fact that the wavelength cannot possibly 
be shorter than the distance between atoms. 

Let this highest frequency be w D . 

In a solid with N identical atoms, there are 3iV energy states, or 
normal modes. This is when we count every possible frequency, 
and every possible "direction of vibration" for each frequency. 

This can be proven mathematically, but we shall skip that. We 
can find u D by integrating the density of states: 

fVD , ., fWfl 3Vuj 2 duj 

/ g{u)dw = = 3N 

■10 JO 2tv ! v j ' 
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The Debye Frequency 



"his equation 



we 3Vuj 2 duj 



, g{w)<h, = = 3N 

JO JO 2n z v :> 

can then be solved to find the Debye frequency. Integrating 

gives 



D 



3N 



Solving for oj d gives the Debye frequency: 

/ 6Ar7r 2„3\ X /3 



v v 



Debye's theory for Heat capacity of solid 
"he phonon distribution is given by 



n{uj)du> 



exp(^/fc B T) - 1 
We can multiply by hw to get the energy at this frequency 
interval: 

riLuniujjdio = ; . 

exp(huj/k B T) - 1 



U = I hiun(uj)dio = / 

Jo Jo 



"he total energy is then 

f u D hwg(u))duj 

exp(!iuj/k B T) - 1' 
We have found earlier that the density of states is: 

r ^ o v " 2 
g{uj) = 3 x 

Substituting into the integral for the total energy above 
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The phonon distribution 



As we have shown earlier, the distribution of phonons is given 
by 

, , , g(s)ds 

n{e)de = ; — , 

exp( £ /fc B T) - 1 

We shall express this in terms of frequency u as well. The 

phonon energy e is given by 

e = fko. 

This comes from assuming that the atoms vibrate in a 
harmonic potential, so that the energy is quantised in steps of 
hu). 

In terms of the frequency w, the phonon distribution is then 
given by 

g(ui)duj 



n(u))du 



exp(hu/k B T) - 1" 
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Debye's theory for Heat capacity of solid 



We get: 

_ 3VH r^D Lu 3 du> 

~ 2tt 2 v 3 JO exp(hu)/k B T) - 1 
At low temperature, it can be shown that this is proportional to 
T 4 . This would imply that the heat capacity C is proportional 
to T 3 . (This has been mentioned in the topic on the specific 
heat of electrons in metal.) 

The T 4 behaviour is derived by making the substitution 

x = huj/k B T. 

As T gets small, x becomes large. So the upper limit is 
substituted with oo. Then using the formula on the last page of 
these notes, it can be integrated to give the T 4 behaviour. 

At high temperatures, the heat capacity tends to 3Nk B . 
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Is verified at high temperature 



In 1912, Debye measured the high temperature behaviour of 
copper. This picture shows a sketch of the molar heat capacity 
of copper that he measured. 



c 

3R 



1.0 



0.5 



_i_ 



_i_ 



100 200 300 

(Deg K) T 



He showed that the heat capacity did tend to 3Nk B , as he had 
predicted. 
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Photons 
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And verified at low temperature 



In 1953, at Purdue University, Keesom and Pearlman measured 
the low temperature behaviour of potassium chloride. This 
picture is a sketch of the results. 

C/T versus T 2 for KCI 



t j 

0 

CD 




As predicted by Debye, the heat capacity was indeed 
proportional to T 3 . 
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Radiation from a very hot object 



An object at a certain temperature would emit and absorb 
electromagnetic radiation. The higher the temperature, the 
higher the average frequency of the radiation emitted. 

When an object gets very hot, this frequency can go as high as 
visible light. For example, the temperature of a volcano lava 
flow can be estimated by observing its color. The result agrees 
well with the measured temperatures of lava flows at about 
1,000 to 1,200 deg C. 




http : //en . wikipedia. org/wiki/Black_body 



The amount of radiation emitted depends on the nature of the 
object - its colour, whether it is smooth or rough, etc. 
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An ideal black body 



We are particularly interested in an ideal black body, one that 
absorbs all the radiation that falls on it. The reason is that the 
radiation emitted by such a body does not depend on its nature 
at all. 

No real material can do this. Soot is about the best, absorbing 
all but 3%. 

In 1859 Gustav KirchhofF had a good idea: "a small hole in the 
side of a large box is an excellent absorber, since any radiation 
that goes through the hole bounces around inside, a lot getting 
absorbed on each bounce, and has little chance of ever getting 
out again." 

(http : / / galileo .phys . Virginia . edu/ classes/252/black_body_radiat ion. html) 

The energy of the absorbed radiation reaches equilibrium 
among different frequencies in the cavity. This results in a 
characteristic spectrum that would be emitted again through 
the hole. 
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Photons in a box 



Like phonons, the photon number is not fixed. We get few 
photons when the box is cold, and more photons when it is hot. 

Unlike phonons, there is no upper limit to the frequency (no 
Debye frequency). The wavelength in a solid cannot be shorter 
than the distance between atoms. The electromagnetic wave in 
a box has no such limit. 

Finally, phonons in a solid can have 3 polarisations: 2 
transverse and 1 longitudinal (like sound). Photons can only 
have 2: both transverse (e.g. light). 
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Photons in a box 



We can derive a formula for the spectrum of black body 
radiation by considering the energies of a free photon gas in a 
box (following Kirchhoff's idea). 

Photons arise from oscillations in electromagnetic fields, which 
is also governed by a wave equation. 

The same ideas for phonons can be used, leading to the same 
density of states: 

g{ -^> = ^~2~3' 

where the speed v is now the speed of light c. This is the 
formula before considering polarisations. 
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Counting photon states 



Recall that for phonons, the density of states is give by: 

VLJ 2 dLU 

g{uj)du} = 3 x 

where the factor of 3 comes from the 3 polarisations of a 
phonon, 2 transverse and 1 longitudinal. 

Since photons only have 2 polarisations, both transverse, the 3 
should be replaced by a 2: 

, . , „ VoJ 2 duj 
g(u)du} = 2 x 

2.n c 

where the sound speed v is replaced by light speed c. 
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Distribution of photon number 



Like phonons, each energy state can be occupied by any 
number of photons. 

So it obeys Bose-Einstein statistics. 



The number density is therefore given by the same formula: 

n{e)de = ; — , 

exp(e/k B T) - 1 

So the number of photons in a given frequency interval is 

Vlj dui doj 

n{Lo)dto = 2 x 



2tt 2 c 3 ex<p(hcu/k B T) - 1 
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Example 



The cosmic microwave background is an example of blackbody 
radiation. The spectrum has been measured. Planck's 
radiation formula is then fitted to the spectrum. 



Wavelength [mm] 

1 0.67 



FIRAS data with 400<J errorbars 
2.725 K Blackbody 




5 14 15 

V [/cm] 



Using a temperature value of 2.74 K in the formula was found 
to fit the measurement well. This tells us that deep space has 
a temperature of 2.74 K. 
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Distribution of photon energy 



The energy of a photon is huj. The energy in the frequency 
interval dw is then 

Vuj 2 duj duj 

hion(u))duj = 2 x nuj x =-=- x 

2n 2 c 3 eXQ(huj/k B T) - 1 

The energy is given by: 

Vhw 3 1 
u(lj) = nuin{Lu) = — =-^= ; ■ . 

This is also formula for the spectrum, and is called Planck's 
law for black body radiation. Integrating gives the total energy: 

2t/i,4 



7T Vk'n . 

U = — ^4t 4 

15h 3 c 3 

It can be shown that this leads to Stefan's law of radiation: 

V = aT 4 

for radiation emitted by an object at a temperature T. 
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Example 

In the same way, Planck's formula may be fitted to the sun's 
spectrum. 



Solar Radiation Spectrum 



/ Sunlight at Top of the Atmosphere 
5250'C Blackbody Spectri 




Radiation at Sea Level 



250 500 750 1000 1250 1500 1750 2000 2250 2500 

Wavelength (nm) 
(http : //en . wikipedia . org/wiki/Sunlight ) 



A 5800 K temperature gives a reasonable fit. However, part of 
the fitted curve deviates from measured data. This could be 
due to absorption by atmosphere, or emission from other 
sources. 
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Example 



The radiation spectrum of a higher temperature object peak at 
higher frequency. 




2.73 K Universe frequency 

(http : //www . launc . tased . edu . au/online/sciences/physics/blackbodyl .html) 

For example, human emit radiation with the peak in the 
infrared wavelength 9.3 /im. A hot plate at 400 °C has a peak 
at 4.3 fj,m, with a bit in visible red. That is why it looks red 
hot. 
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Example 1 

(i) Write down the integral expression for the total energy of 
phonons in a solid. Write down the formula for Debye 
frequency. Define all symbols. 

(ii) Make the substitution x = huj/k B T where w is the angular 
frequency and T is the temperature. Find the expression for 
the integral for at the low temperature limit. 

(iii) Integrate to find the total energy. Show that phonon heat 
capacity is proportional to T 3 at the low temperatures. 
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Worked Examples 
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Solutions 

(i) Expression for total energy is 

r u D hujg(ijj)ckjj 



U 



Jo exp(hu}/k B T) - 1 
where u> D is the Debye frequency, w the angular frequency, g(uj) 
the density of states and T the temperature. 

Formula for Debye frequency is 

where iV is total number of atoms, v is speed of sound and V is 
volume of solid. 



(ii) Rearrange: 

k B Tx 

UJ = 

h 



31 



For the upper limit, lj d , of the integral, let the corresponding x 
be x£). So 

hior) 

x d = TT^p- 

k B T 

Density of states of phonons is: 

( \a 3Vuj2 j 



Transforming to x variable 

2 7r 2 u 3 \ 



g(x)dx = , (— - ^ B dx. 



Regrouping terms: 

, , 3Vk%T 3 2 

q{x)dx = _ %, - X dx. 

2n 2 v 3 h 3 

Combining the above into the integral: 

f x D q(x)dx 
U = I k B Tx x 



Jo exp(x) — 1 
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Substitute into the low temperature limit above, we finally get 

3Vk%T* 7T 4 

U 



2-K 2 v 3 h 3 15' 

Differentiating with respect to T, we get the heat capacity: 

„ 2Vk%n 4 , 
C = R 2 3*3 X T ' 

which is proportional to T 3 . 
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Substitute the expression for g(x)dx: 

fx D 3Vk 3 T 3 x 2 dx 
U = / k B Tx x — ^4 q x , 

Factor terms that do not depend on x: 

3Vk%T 4 r*D x 3 dx 

U = — x / • 

2n z v :i h :i Jo e x — 1 

If T goes to zero, 

fc B T 

goes to infinity. So at the low temperature limit, 

3Vk%T A foo x 3 dx 

U = — =4"7f x / . 

2n 2 v 3 h 3 Jo e x - 1 



(iii) To find the total energy, we use the result 

--3 4 

X J 7T 

-dx 



Jo e x - 1 15 
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Example 2 

(i) Write down the integral expression for the total energy of 
photons in a box. 

(ii) Make the substitution x = tiw/k B T where oj is the angular 
frequency and T is the temperature. Integrate to find the total 
energy. 

(iii) Stefan's law gives the power per unit area coming out in all 
directions through a hole on the box ('per unit area' referring 
to area of the hole). Assuming that radiation is only emitted in 
the direction normal to the area, find an approximate formula 
for Stefan's law. Why is this approximate? 
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Solution 



(i) The expression for total energy is 

r°° Vhuj 3 1 

U = / ^-^ ; ■ r du!. 

Jo tt 2 c 3 exp(hu/k B T) - 1 
where V is volume and T is temperature. 

(ii) Rearranging, 

k B T 

LO = X. 

h 

Substituting into the integral: 



h n 2 c 3 \ h J exp(a;) — l h 
The limits remain the same since x and uj are proportional 
Regrouping, 

Vk%T A roo x 3 
U = ^4 „ x / da;. 



7r 2 c 3 fi 3 Jo e 31 — 1 
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The energy density is energy / volume = U/V. So we have 

_ Uc 

Substituting the above result for U, we get 

_^ Vk%T A 
V ~ V X ir 2 c 3 h 3 x 15' 
Simpifying, we have an approximate Stefan's law: 

7T 4 k 4 
15TT 2 C 2 h 3 

This is approximate because radiation can really come out in all 
directions. We have not taken this into account. 
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Using this result: 

dx = — . 
15 



r°° x~" 
Jo e x - 1 



we find 



Vk%T 4 ^ 
U = b „ x — . 
TY z c~ > n~ > 15 



(iii) Recall the formula for water flowing through a pipe: 

mass per second = density x area x velocity. 

Instead of mass, we have energy here. So the formula becomes 

power = energy density x area x speed of light. 

For Stefan's law, we want the intensity r\\ 
power 

= energy density x speed of light. 

area 
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Learning Outcome: You should be able to 



State and use the phonon energy distribution. Derive the 
density of states for phonons. Explain and derive the Debye 
frequency. 

Derive the low and high temperature limits of the phonon heat 
capacity. 

Explain blackbody radiation as a free photon gas. State and 
use the photon energy distribution. Derive the density of states 
for photons. 

Derive Stefan's law for blackbody radiation. 
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Useful Integrals 
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Useful Integrals 



For Debye model, to find heat capacity at low temperature; 
and for black body radiation, to find Stefan's law: 



rOC x ° 

JO e x - 1 



dx = — . 
15 



For liquid helium-4, later in the course:- 

to find condensation temperature of Bose Einstein condensate: 



roc X '-i'- 
Jo e x - 1 

to the find heat capacity: 



dx = 2.315 



roc x^' 2 

/ — dx= 1.783 

Jo 



10 e x - 1 
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